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ABSTRACT 

Context. Convection in red giant stars excites resonant acoustic waves whose frequencies depend on the sound speed inside the star, 
which in turn depends on the properties of the stellar interior. Therefore, asteroseismology is the most robust available method for 
probing the internal structure of red giant stars. 

Aims. Solar-like oscillations in the red giant HR 7349 are investigated. 

Methods. Our study is based on a time series of 380 760 photometric measurements spread over 5 months obtained with the CoRoT 
satellite. Mode parameters were estimated using maximum likelihood estimation of the power spectrum. 

Results. The power spectrum of the high-precision time series clearly exhibits several identifiable peaks between 19 and 40 ^Hz 
showing regularity with a mean large and small spacing of Av = 3.47 ± 0.12yuHz and (5vo2 = 0.65 ± 0. 10/jHz. Nineteen individual 
modes are identified with amplitudes in the range from 35 to 1 15 ppm. The mode damping time is estimated to be days. 

Key words. Stars: individual: HR 7349 - Stars: oscillations ~ Stars: variables: general - Stars: interiors 



1. Introduction 

' The analysis of the oscillation spectrum provides an unrivaled 
' method for probing the stellar internal structure because the fre- 
quencies of these oscillations depend on the sound speed in- 
side the star, which in turn depends on the density, temperature, 
gas motion, and other properties of the stellar interior. High- 
precision spectrographs have acquired data yielding to a rapidly 
] growing list of solar-like oscillation detections in main-sequence 
. and giant stars (see e.g.. Bedding & Kieldsen |200 7| Carrier et al. 
' I2008I I. In a few years, we have moved from ambiguous detec- 
' tions to firm measurements. Among these, only a few are re- 
lated to red giants, e.g., ^ Hya, Frandsen et al. (2002), e Oph, De 
Ridder et al. (2006 ), and rj Ser, Barban et al. (2004). The reason 
is that longer and almost uninterrupted time series are needed to 
characterize the oscillations in red giants, because of longer os- 
cillation periods than main-sequence stars, and long observing 
runs are difficult to obtain using high-accuracy spectrographs. 

The CoRoT (COnvection Rotation and planetary Transits) 
satellite (Baglin 2006) is perfect for this purpose because it 
can provide these data for a large number of stars simulta- 
neously. The CoRoT satellite continuously collects white-light 
high-precision photometric observations for 10 bright stars in 

* The CoRoT space mission has been developed and is operated 
by CNES, with the contribution of Austria, Belgium, Brazil, ESA, 
Germany and Spain. 



the so-called seismofield, as well as 3 -color photometry for thou- 
sands of relatively faint stars in the so-called exofield. The pri- 
mary motivation for acquiring this second set of data is to de- 
tect planetary transits, but the data are also well suited to as- 
teroseismic investigations. De Ridder et al. C2009') unambigu- 
ously detected long-lifetime non-radial oscillations in red giant 
stars in the exofield data of CoRoT, which is an important break- 
through for asteroseismology. Indeed, observations from either 
the ground or other satellites have been unable to confirm the 
existence of non-radial modes and determine a clear value of the 
mode lifetime. Hekker et al. (I2009I I presented a more detailed 
classification of the red giants observed by CoRoT. 



HR 7349 (HD 1 8 1 907) is a bright equatorial G8 giant star (V 
= 5.82) that is an excellent target for asteroseismology. This star 
was selected as a secondary target during the first long run of 
the CoRoT mission. In this paper, we thus report on photomet- 
ric CoRoT observations of HR 7349 resulting in the detection 
and identification of p-mode oscillations. The non-asteroseismic 
observations are presented in Sect. 2, the CoRoT data and fre- 
quency analysis in Sects. 3 and 4, and the conclusions are given 
in Sect. 5. 
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2. Fundamental parameters 

2.1. Effective temperature and chemical composition 

We used the line analysis code MOOG, Kurucz models, and a 
high-resolution FEROS spectrum obtained in June 2007 to cany 
out an LTE abundance study of HR 7349. The effective temper- 
ature and surface gravity were estimated from the excitation and 
ionization equilibrium of a set of iron lines taken from the line 
list of Hekker & Melendez (2007). We obtain Teff =4790±80 K 
and [Fe/H]=-0.08±0.10 dex, while the abundance pattern of the 
other elements with respect to Fe is solar within the errors. The 
full results of the abundance analysis will be reported elsewhere 
(Morel et al., in preparation). We also determined a photomet- 
ric temperature given by the relation in Alonso et al. dl 9991 1 us- 
ing the dereddened color index (B-V) (see Sect. I2.2l i and found 
4704+llOK, which agrees with the spectroscopic value. We fi- 
nally adopt a weighted-mean temperature of 4760+65 K. 

2.2. Luminosity 

Even for such a bright star, the interstellar extinction in the direc- 
tion of the Galactic center is not negligible. From the Hipparcos 
parallax 11 - 9.64 ± 0.34 mas (van Leeuwen 120071 1 and the 
value of (B-V)= 1.093 in the Hipparcos catalog, an absorption 
of Ay = 0.1 85 mag is derived for the region of the star (Arenou 
et al. 1199 2). which corresponds to an Eb-v = 0.052. Combining 
the magnitude V - 5.809 + 0.004 (Geneva photometry, Burki 
et al. 120081 ). the Hipparcos parallax, the solar absolute bolomet- 
ric magnitude Mboi,G - 4.746 (Lejeune et al. ,199 8). and the 
mean bolometric corrections EC = -0.40 ± 0.04 mag (BC - 
-0.42 + 0.06 mag according to the calibration of Flower, 119961 
and BC = -0.38 + 0.05 mag by Alonso et al., ,1999i ), we find a 
luminosity for HR 7349 of L = 69 ± 6 Lq. 

2.3. Rotational velocity 

We determined the rotational velocity of the star by means of 
a spectrum taken with the spectrograph Coralie installed on 
the 1.2-m Swiss telescope, ESO La Silla, Chile. According 
to the calibration of Santos et al. (120021) . we determined a 
V sin / = 1 .0 + 1 .0 km s ' . For this small value of the projected ro- 
tational velocity, we do not expect to see any split modes in the 
power spectrum. 

2.4. Large spacing estimation 

An estimation of the mass of HR 7349 may be obtained by 
matching evolutionary tracks to the L-T^ff error box in the HR di- 
agram. However, in the red giant part of the HR diagram, this de- 
termination is not robust at all. Assuming a mass of between 0.8 
and 3 and scaling from the solar case (Kjeldsen & Bedding 
119951 ). a large frequency spacing of 2.8-5.5 //Hz is expected. 



3. CoRoT Observations 

HR 7349 was observed with the CoRoT satellite for 5 consecu- 
tive months. CoRoT was launched on 2006 December 27 from 
Baikonur cosmodrome on a Soyuz Fregat 11- lb launcher. The 
raw photometric data acquired with CoRoT were reduced by the 
CoRoT team. A detailed description of how photometric data 
are extracted for the seismology field was presented in Baglin 
( 120061 ). A summary can be found in Appourchaux et al. (i2008i) . 
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Fig. 1. The total CoRoT light curve (top) and a zoom (bottom) 
of HR 7349. This light curve is detrended with a polynomial fit 
(order 8), which only affects frequencies below 6 yuHz. A period- 
icity of about 8.5 hours can be seen in the zoom, corresponding 
to oscillation modes close to 30yuHz. 

In the seismofield, CoRoT obtains one measurement ev- 
ery 32 seconds. The observations lasted for 156.64 days from 
May 11th to October 15th 2007. The light curve shows near- 
continuous coverage over the 5 months, with only a small num- 
ber of gaps due mainly to the passage of CoRoT across the South 
Atlantic Anomaly. These short gaps were filled by suitable in- 
terpolation (Baglin 120061) . without any influence on the mode 
extraction because it only affects the amplitude of frequencies 
far above the oscillation range of our target (see Fig. |2|. The 
duty cycle for HR 7349 before interpolation was 90 %. For the 
frequency analysis (see Sect. H)), the light curve was detrended 
with a polynomial fit to remove the effect of the aging of the 
CCDs (see Auvergne et al. 2009). This detrending has no con- 
sequence on the amplitude or frequency of oscillation modes, 
since it only affects the power spectrum for frequencies lower 
than 6 fjHz. The light curve shows variations of a timescale of 8- 
9 hours and peak-to-peak amplitudes of l-3mmag (see Fig.[T]). 
This signal is a superposition of tens of smaller modes with sim- 
ilar periods (see Sect.|4]). 

4. Frequency analysis 

4.1. Noise determination 

We computed the power spectra of CoRoT light curves with 
both gaps and interpolated points (see Sect. |3]). The resulting 
power spectra are quasi-identical, the interpolation not affect- 
ing the oscillations but suppressing the aliases (the most impor- 
tant of which lies at 23 yuHz). We thus analyze the interpolated 
time series with negligible alias amplitudes. The time base of the 
observations gives a formal resolution of 0.07 /iHz. The power 
(density) spectrum of the time series, shown in Figs. |3] and |4] 
exhibits a series of peaks between 20 and 40 //Hz, exactly where 
the solar-like oscillations are expected for this star. The power 
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Fig. 2. Power spectra of the original data (grey) and interpolated 
data (black). The range of the oscillation is zoomed in the inset. 
The interpolation drastically reduces the amplitude of aliases, in 
particular the one at 23 //Hz. 
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Fig. 3. Power density spectrum of the photometric time series 
of HR 7349 and a multi-component function (black line) fitted 
to the heavily smoothed power density spectrum. The function is 
the superposition of three power-law components (dashed lines), 
white noise (horizontal dashed line) and a power excess hump 
approximated by a Gaussian function. 
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Fig. 4. Power spectrum of the CoRoT observations of HR 7349. 
Only a polynomial fit was removed from the original data. 



law. For a given component, the power remains approximately 
constant on timescales longer than B,, and drops off for shorter 
timescales. Each power law corresponds to a separate class 
of physical phenomena, occurring on a different characteristic 
timescale, and corresponding to different physical structures on 
the surface of the star. In our case, we fixed the slope to 4, which 
is a typical value for the Sun (Aigrain et al. 120041 Michel et al. 
,2009). Moreover, this value allows us to fit our power density 
spectrum well. 

To model the power density spectrum, we added a white 
noise P„ and a power excess hump produced by the oscillations, 
which was approximated by a Gaussian function 



A ■ 

^)=i;tt(b^^^"^^^ 



-(v„„-v)"/(2<r2) 



(2) 



density spectrum is independent of the observing window: this is 
achieved by multiplying the power by the effective length of the 
observing run (we have to divide by the resolution for equidis- 
tant data), which is calculated to be the reciprocal of the area 
beneath the spectral window in power (Kjeldsen et al. l2005l l. We 
note that to obtain the same normalisation as in Baudin et al. 
(I2005I I. we multiply the power by the effective length of the ob- 
servation divided by four. Typically for such a power spectrum, 
the noise has two components: 

- At high frequencies it is flat, indicative of the Poisson statis- 
tics of photon noise. 

- Towards the lowest frequencies, the power scales inversely 
with frequency, as expected for instrumental instabilities and 
noise of stellar origin like granulation. 

For the Sun, it is common practice to model the background sig- 
nal with power laws to allow accurate measurements of solar os- 
cillation frequencies and amplitudes (Harvey 1985 Andersen et 
al.[i998 , Aigrain et al. 2004). To study this "noise", we compute 
the power density spectrum shown in Fig. [3] and fit a smoothed 
version of this spectrum with a sum of N power laws 



P(y) 



1=1 1=1 



(1) 



where this number N depends on the frequency coverage, v is 
the frequency. A, is the amplitude of the i-th component, B, 
is its characteristic timescale, and C, is the slope of the power 



The number of components N is determined iteratively: we first 
made a single component fit, and additional components were 
then added until they no longer improved the fit. In our case, we 
limited the number of components to three. We note that for this 
star, part of the low frequency variation is caused by the aging 
of the CCD. The timescales for the different noise sources (B,) 
are 3.6 x 10^ 6.9 x 10"^, and 1.1 x 10"* s. The noise at high 
frequencies P„ is only 0.17 ppm^/yuHz and the oscillations are 
centered on 28.2 juHz. 



4.2. Search for a comb-like pattern 

In solar-like stars, p-mode oscillations are expected to pro- 
duce a characteristic comb-like structure in the power spectrum 
with mode frequencies v„ f reasonably well approximated by the 
asymptotic relation (Tassoul 1980): 



v„,f ^ Av(« +- + e)-e(e+ l)Do 



(3) 



Here Do (which equals ^5vo2 if the asymptotic relation holds 
exactly) and e are sensitive to the sound speed near the core 
and in the surface layers, respectively. The quantum numbers 
n and { correspond to the radial order and the angular degree 
of the modes, and Av and 6vo2 are the large and small spac- 
ings. We note that a giant star such as HR 7349 is expected to 
show substantial deviations from its regular comb-like structure 
described above (Christensen-Dalsgaard 2004). This is because 
some mode frequencies, except for { - 0, may be shifted from 
their usual regular spacing by avoided crossings with gravity 
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Fig. 5. Autocorrelation of the slightly smoothed power spectrum. 
The large spacing, as well as separations between { = and (-1 
modes are clearly present. 



modes in the stellar core (also called 'mode bumping') (see e.g., 
Christensen-Dalsgaard et al. 11995' and Fernandes & Monteiro 
I2003I I. We must keep the possibility of these mixed modes in 
mind when attempting to identify oscillation modes in the power 
spectrum. Moreover, the ratio of lifetime to oscillation periods is 
usually far smaller for red giants (because of their longer pe- 
riods) than for solar-like stars, which can complicate the mode 
detection (see e.g., Stello et al. 120061 TaiTant et al. 2008). 

The first step is to measure the large spacing that should ap- 
pear at least between radial modes. The power spectrum is auto- 
correlated to search for periodicity. Each peak of this autocorre- 
lation (see Fig.|5]l corresponds to a structure present in the power 
spectrum. One of the three strong groups of peaks at about 1 .7, 
3.5, and 5.2yL(Hz should correspond to the large spacing. By vi- 
sually inspecting the power spectrum, the value of about 3.5 jjHz 
is adopted as the large separation, the two others are spacings be- 
tween the { = and f-l modes. This large separation value is 
in good agreement with the scaled value from the solar case (see 
Sect. 12311. 

4.3. Extraction of mode parameters 

The power spectra in Figs. |6]and |7] clearly exhibit a regular- 
ity that allows us to identify { = to {-2 modes. Each mode 
consists of several peaks, which is the clear signature of a fi- 
nite lifetime shorter than the observing time span. In order to 
determine the mode frequencies, as well as amplitude and life- 
time of the modes, we fitted the power spectrum using a maxi- 
mum likelihood estimation (MLE) method. MLE has been ap- 
plied widely in the helioseismic community (e.g., see Schou 
[19921 Appourchaux et al. 1998 , and Chaplin et al. 2006). Our 
program uses the IDL routines developed by T. Appourchaux 
( I19981 I. This method has already been used with success for the 
red giant e Oph (Barban et al. 2007). 

The modelled power spectrum for a series of M oscillation 
modes, P(vjt), is 



M 



H„ 



1 



(4) 



where H„ is the height of the Lorentzian profile, v„ is the os- 
cillation mode frequency, F,, is the mode line-width, and B is 
the background noise. The fit was performed on the non-over- 




I'Vcqucncy | /iWv. 

Fig. 6. Power spectra of the whole dataset (black) and of the two 
half-long subsets (grey). The frequencies of the oscillation peaks 
change from one subset to the other, which is a sign of a finite 
lifetime. The initial guess for the identification of modes are in- 
dicated by shaded regions: these regions are regularly spaced in 
agreement with the large separation deduced from the autocor- 
relation. Every second region has a structure more simple and 
narrow, corresponding to our identification of ^ = 1 modes, the 
others are more complex and wide and correspond to ^ = and 2 
modes. 



sampled power spectrum to minimize the interdependency of the 
points. The quantity that was minimized is 



L = ^ lnP(y,) + 



PobsiVk) 
PiVk) 



(5) 



where K is the number of bins, i.e., the number of Fourier fre- 
quencies. 

Only a few peaks belong to a given mode, indicating a life- 
time longer than 10 days. It is thus also difficult to derive the 
correct Lorentzian shape for each mode. Therefore some param- 
eters were fixed to avoid incorrect parameter determinations: 

- The noise was determined independently of the MLE method 
(see Sect.l4Tb. 

- Since no difference between the shapes of modes of dif- 
ferent degree t is detected and the value of vsin/ is ex- 
tremely small, resulting in an expected rotational splitting 
value smaller than 0.03 yuHz, we assumed that it is zero. 

- When fitting all modes without fixing the width of the 
Lorentzian envelope, we clearly saw that the fit was not ro- 
bust enough to provide an accurate determination of all mode 
parameters. The method was thus first to find a mean value 
for the Lorentzian width and to fix this mean value for all 
modes. The mean width of the Lorentzian is obtained by in- 
dividually fitting all modes and by taking their mean, reject- 
ing the values that obviously correspond to a poor fit. The 
determined mean is 0.25 ± 0.06 yuHz, which corresponds to a 
mode lifetime of 14.7+^-^ d. The use of sub-series show the 
stochastic nature of the mode excitation caused by the dif- 
ferent fine structure of peaks in each spectrum, and give an 
indication of their width. We note that we also checked that 
no significant width difference was found for modes with dif- 
ferent degrees d (by comparing their mode-width means). 
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Fig. 7. Top: Lorentzian fit (thick red line) of the observed power spectrum assuming the same lifetime for all the oscillation modes. 
Bottom: Ratio of the observed spectrum over the Lorentzian fit (in amplitude). As expected, this ratio does not show any correlation 
with the fit. 



All modes were then fitted with fixed parameters as deduced 
above. We note that the initial guess values are indicated by 
shaded regions in Fig. |6] These regions are already in good 
agreement with a regular spacing between modes. Every sec- 
ond region has a structure that is more simple and narrow, cor- 
responding to our identification of f = 1 modes, the others are 
more complex and wide. In this last case, we needed to fit two 
modes per region to reproduce the power spectrum (which cor- 
respond to { = and 2 modes). At 38 yuHz, the signal-to-noise 
ratio is far smaller and it becomes more difficult to differentiate 
between modes: for the sake of homogeneity, we also decided 
to fit two modes in this region. The frequencies of these last 
two modes will however be more uncertain. The result is listed 
in Table [1] The formal uncertainties associated with MLE are 
well understood, as explained by Libbrecht (I1992i l and Toutain 
& Appourchaux (.1994) . The echelle diagram with the nineteen 
identified modes is shown in Fig. [8] At higher and lower fre- 
quency (above 40 //Hz and below 19 //Hz), the amplitude of the 
modes is either too small or the noise too high to unambiguously 
identify additional modes. The ratio of the observed power spec- 
trum to the fit is shown in Fig. [T] it appears to be pure noise 
and has a mean value of 1 . Moreover, no correlation was found 
between the amplitude of this ratio and the amplitude of the fit. 
We note that the non-radial modes are as well aligned as radial 
or non-mixed modes. However, the separation between {-I and 
{ =0 modes is not fully compatible with the asymptotic relation. 
The {=1 modes are indeed too far to the right in the echelle dia- 
gram. 

The large and small separations are shown in 
Fig. |9l The mean large separation has a value of 
Av = 3.47 ±0.1 2 //Hz, and the values for different degrees 
{, Ave, are: Av,, = 3.45 ± 0.12//Hz, Avi = 3.46 + 0.07 //Hz, and 
Av2 = 3.50 ± 0.19//Hz. We can identify a small oscillation 
of the large spacing that varies with frequency, which is a 
clear signature of the second helium ionization zone (see e.g., 
Monteiro & Thompson i 19981) . The small separation has a mean 
value of 5v()2 = 0.65 ± 0.10 //Hz and seems to decrease with 
frequency. The small value of the frequency difference between 
these modes make their frequency determination more uncertain 
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Fig. 8. Echelle diagram of identified modes with a large separa- 
tion of Avo = 3.45 //Hz. The modes { = {•),( = 1 (a), and ^ = 2 
(o) follow ridges, but (-1 modes are situated too far to the right 
to follow the asymptotic relation. 

than that of £- \ modes, which are not affected by neighbouring 
modes. 



4.4. Oscillation amplitudes 

The fit of the Lorentzian profiles to the power spectrum infers 
the height of all oscillation modes. Since the modes are resolved 
and because of the normalization of the power spectrum, the rms 
amplitude is measured to be (see Baudin et al. 120051 ) 



(6) 



where H and F are the height and width (FWHM) of the 
Lorentzian function, respectively, in the power density spectrum. 
The amplitudes are in the range of 35 - 1 15 ppm (see Fig. [TOl i. 
The error bars are derived from the Hessian matrix and the corre- 
lation between the height and width has been taken into account 
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Table 1. Frequency and amplitude of identified modes. The un- 
certainties in the last given digit of the frequencies are noted in 
parentheses. For the S/N estimates (see Barban et al. I2007I I. the 
signal is taken to be the height of the fitted Lorentzian profile and 
the noise is determined according to the procedure described in 
Sect.Q] 
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Fig. 9. Top: Small spacing versus frequency. Bottom: Large 
spacing versus frequency. The variations of the large separation 
with frequency show a clear oscillation. The symbols used are 
the same as for Fig. [8] 



(Toutain & Appourchaux 1 1 9941 Appourchaux, private commu- 
nication). We note that the amplitudes of the l- \ modes are the 
highest. 

As noticed by Kjeldsen et al. (|2005l). measurements made 
on different stars with different instruments using different tech- 
niques, in different spectral lines or bandpasses, have different 



sensitivity to the oscillations. It is thus important to derive a 
bolometric amplitude that is independent of the instrument used. 
We computed the maximum bolometric amplitude of the /" = 
modes, because their visibility coefficients do not depend on the 
inclination of the star. According to Michel et al. (I2009I I. who de- 
rived the CoRoT response for radial modes, radial-mode ampli- 
tudes of HR 7349 must be divided by 1 . 1 6 to obtain the bolomet- 
ric amplitudes. We find that Aboi/=o,max = 81 ppm, which corre- 
sponds to 32 times the solar value (Michel et al. 120091 ). The scal- 
ing laws for both the large separation and the frequency of the 
maximum amplitude (Kjeldsen & Bedding, 11 9951 ), coupled with 
the non-asteroseismic constraints, infer a mass for HR 7349 of 
about 1 .2 Mq The derived amplitude is in good agreement with a 
scaling function (L/M)V -y/Tdf/T^efi^o, with s close to 0.8, which 
is in-between the values given by Samadi et al. (120071 s=0.7) 
and Kjeldsen & Bedding (fT995l s=l). 



5. Conclusion 

The red giant star HR 7349 has been observed for about 156 
days by the CoRoT satellite. These observations have yielded a 
clear detection of p-mode oscillations. As already mentioned by 
De Ridder et al. (2009), non-radial modes are observed in red 
giants. Nineteen identifiable modes of degree l-^\al-2 ap- 
pear in the power spectrum between 1 8 and 42 juHz with average 
large and small spacings of 3.47 and 0.65 yuHz, respectively, and 
a maximum bolometric amplitude of 8 1 ppm. We note that the 
amplitude of ^ = 1 modes is even larger than that of the radial 
modes. 

All modes of the same degree are aligned in the echelle di- 
agram, which is a sign of modes that follow the asymptotic re- 
lation. However, the separation between (,—\ and I =0 modes is 
not fully compatible with this asymptotic relation. All frequency 
patterns are theoretically expected for red giants (Dupret et al. 
120091 ), from very complex to regular: our observations corre- 
spond to a red giant for which the radiative damping of non- 
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radial modes is large and only radial and non-radial modes com- 
pletely trapped in the envelope can be observed. 

By fitting Lorentzian profiles to the power spectrum, it has 
also been possible to unambiguously derive, for the first time for 
a red giant, a mean line-width of 0.25 fiHz corresponding to a 
mode lifetime of 14.7 days. This lifetime is in agreement with 
the scaling law T^^ suggested by Chaplin et al. I2009i although 
is a little too long. This relation however has yet to be verified for 
a larger number of red giants with different physical properties. 
The theoretical study of this red giant, including asteroseismic 
and non-asteroseismic constraints, will be the subject of a second 
paper. 
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